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Oh 1 Abstract 

O ; 

^ \ Under certain circumstances, the group velocity in a homogenized composite medium 

(HCM) can exceed the group velocity in its component material phases. We explore 
this phenomenon for a uniaxial dielectric HCM comprising isotropic component ma- 
terial phases distributed as oriented spheroidal particles. The theoretical approach is 

Qh' based upon the Bruggeman homogenization formalism. Enhancement in group velocity 

in the HCM with respect to the component material phases is shown to be sensitively 
dependent upon the shape of the component spheroids and their alignment relative to 
^ \ the direction of propagation. 
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' 1 Introduction 

£* 

The process of homogenization involves the combination of two (or more) component material 
^ ■ phases to produce a single, effectively homogeneous, composite medium [1, 2, 3]. Typically, 
the constitutive properties of the component material phases are relatively simple as com- 
pared with those of the homogenized composite medium (HCM). Through homogenization, 
novel and potentially useful material properties may be realized [4, 5]. Many examples of 
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material properties being extended — or indeed entirely new material properties being real- 
ized — as a result of homogenization can be found within the regimes of linear and nonlinear 
electromagnetics [6]. 

An interesting result concerns the electromagnetic group velocity in HCMs. Under cer- 
tain circumstances, the group velocity in an HCM can exceed the group velocities in its 
component material phases. This issue has been investigated for isotropic dielectric com- 
posite mediums using the Maxwell Garnett [7, 8] and the Bruggeman [9] homogenization 
formalisms. In these studies, an enhancement in group velocity is demonstrated through 
homogenizing two component material phases, one of which is characterized by a relatively 
large permittivity and relatively small frequency-dispersive term as compared with the other 
component material phase. 

Enhancement of group velocity in a laminate composite medium has been considered by 
using a volume-weighted sum to estimate the HCM permittivity [7, 8]. The directional prop- 
erties of group- velocity enhancement are further explored in this communication. Specifi- 
cally, we consider a uniaxial dielectric HCM which develops from the homogenization of a 
random assembly of oriented spheroidal particles. The component material phases are them- 
selves electromagnetically isotropic. Our theoretical analysis is founded upon the Bruggeman 
homogenization formalism [10]. 

2 Homogenization 

Let us consider the homogenization of a composite medium containing two component ma- 
terial phases, labelled as a and b. Both component material phases are taken to be isotropic 
dielectric mediums: e a and e b denote the permittivity scalars of phases a and b, respectively. 
In order to focus in particular upon the phenomenon of enhancement of group velocity, 
without being distracted by the complications arising from dielectric loss, the component 
material phases are assumed to be nondissipative; i.e., e a ' b e M. The component material 
phases are envisioned as random distributions of identically oriented, spheroidal particles. 
The spheroidal shape — which is taken to be the same for all particles of phases a and b — is 
parameterized via the shape dyadic 



where £ is the identity 3x3 dyadic and the unit vector c is parallel to the spheroid's axis of 
rotational symmetry. The spheroid's surface is described by the vector 



with f being the radial unit vector from the spheroid's centroid and specified by the spherical 
polar coordinates and <fi. The linear dimensions of the spheroid, as determined by the 
parameter 77, are assumed to be small relative to the electromagnetic wavelength (s). 
The permittivity dyadic of the resulting HCM, 



g=U ± L + (U\\-U±) cc 



(1) 



r 8 {e,4>) = 7 l U't{e,4>), 



(2) 




Br 




(3) 
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is estimated using the Bruggeman homogenization formalism as the solution of the equation 

f a g a + f b a b = g, (4) 

where /„ and fb = 1 — f a denote the respective volume fractions of the material component 
phases a and b [10]. The polarizability dyadics in (4) are defined as 

i = {e'L - g Br ) • [L + R' {e'L ~ g Br ) ]~\ (t = a,b), (5) 
wherein the depolarization dyadic is given by the surface integral [11, 12, 13] 

D = — [ d(f) I dOsmel , 1 , | U' 1 -ff -U~ l . (6) 

The depolarization dyadic may be expressed as 

R = D ± L+{D ll -D ± )cc, (7) 

where 

D\\ = 71^(7), (8) 

D± = 4^(7), (9) 
e _i_ 

The terms Tu and T± herein are functions of the real- valued parameter 



U 2 ± ef r 
Uh 1 



7 = 772i; (io) 



they have the representations 



rw n 1 f , , f m cos 2 0sin 3 fl , , 

r,, 7 = — d6 d6 — =— f = 5— r, (11) 

l|V ; 4tt J y J cos 2 9 + sin 2 (7 cos 2 + sin 2 0) V ; 

, . 1 f 2n T sin 2 0sin 3 # . , 

r_L 7 = — d<b d6 — =— £ 5-^-. (12) 

w; 4tt 7 r 7 cos 2 + sin 2 (7 cos 2 + sin 2 0) V ; 



3 



The surface integrals (11) and (12) may be evaluated as 

1 



r,i(7) 



sinh 



-1 / 1-7 

7 



(l- 7 )3 1-7 
1 sec" 1 ^/7 



I 7-1 ( 7 -l) 



1 / 1 



7 sinh 



-1 /l- 7 

7 



r±(7) = < 



mi -7 (i_7) 

1 / 7 sec" 1 y/7 1 
( 7 -l)I 7-l y 



for < 7 < 1 
for 7 > 1 
for < 7 < 1 

for 7 > 1 



(13) 



(14) 



We exclude the cases of 

• the isotropic HCM with 7 = 1, and 

• the anomalous hyperbolic HCM with 7 < [14] 
from consideration. 

The dyadic Bruggeman equation (4) provides the two nonlinear scalar equations 



-/« + 



e b - ef 



1 + D\\ (e a - ef r ) ° 1 + D\\ (e 6 - ef r ^j 



h = o, 



c a ,Br 
t t 1 



1 + D± (e a - e 



Br 



■/« + 



€ — € 



Br 



1 + D ± (e b - ef) 



■fb = 0, 



(15) 



(16) 



coupled via D± t \\, which can be solved straightforwardly for ef and ef using standard 
numerical techniques. 

3 Group velocity 

Let us consider a wavepacket which is a superposition of planewaves with phasors 



K{r) = E_ exp (ik • r) 
H_{r_) = H_q exp (ik • r) 



(17) 



The group velocity v g of the wavepacket is conventionally defined in terms of the gradient 
of the angular frequency u with respect to k [15]; i.e., 



v = V fc u 



Sl luj=Lu(k a vg) 
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(18) 



where k avg denotes the average wavenumber of the wavepacket. Herein we adopt the compact 
notation 

V - (— — JL\ (iQ) 
- \ dk x ' dk y ' dk z J 

for the gradient operator with respect to k, where (k x ,k y ,k z ) is the representation of k in 
terms of its Cartesian components. 

In order to calculate the group velocity in the uniaxial dielectric HCM (3), denoted as 
v_g r , we exploit the corresponding planewave dispersion relation as follows. The combination 
of (3) with the source-free Maxwell curl postulates 



V x E(r) = iuoB(r) 

V x H(r) = -iujD(r) 



(20) 



delivers the vector Helmholtz equation 

[(Vxl) • (Vxlj-^f] 'E = 0, (21) 

with yU being the permeability of free space. The requirement that (21) provide nonzero 
solutions for the planewave phasors (17) yields the dispersion relation 

W(k,u)=0, (22) 

wherein the scalar function W is defined as 

W(k, u) = (k - k - ef r fi uj 2 ) (k • f r • k - ef r ef r fi u 2 ) . (23) 

The dispersion relation (22) admits two wavevector solutions: the ordinary wavevector k 
and the extraordinary wavector k ex , satisfying 



(24) 



We note that the magnitude of the ordinary wavevector is direction-independent, and the 
ordinary and extraordinary wavevectors coincide when k ex is directed along c. 

By taking the gradient of the dispersion relation (22) with respect to k, we find 

dW 

V k W+—\7 k iu = 0. (25) 
ouj 

Hence, the HCM group velocity (18) may be expressed as 

1 



9 dW/du ~ 



(26) 

u}=u)(k avg ) 
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The partial derivative terms involving W are found to be 
V k W = 2[(k-g Br 'k-uo 2 fi e B ; r ef r )k+(k'k-uo 2 fi ef r )g Br -k], 



(27) 



dW 

duo 



{k m k — cJ 2 /i ef r ) ^ k • -| — • k — jj, uo 



Br 



2e Br e Br + uo 



duo ± duo 



W (2ef + uo^p) (k'g Br • k — ef ef /W) , 



with 



de Br de Br -j. /def de 



Br 



do; do; 



/ 



da; duo 



cc. 



(28) 



(29) 



By virtue of (24), we see that the ordinary and the extraordinary group velocities are given 
by 

2 



.Br I 



(30) 



and 



Br I 



o;/i 



2ef ef r + o> 



<fef Br , «r dc 



Br 



duo 



duo 



de Br 



(31) 



respectively. 

In order to find the derivatives of ef 1 " and e Br needed to evaluate the group velocities 
(30) and (31), we have to exploit the Bruggeman equations (15) and (16). As a precursor, 
let us first note the derivatives of the depolarization dyadic components 

dD n def de Br 

'- ; h«i2- 



duo 
dD ± 
duo 



= OL 2 \- 



duo 
def 

duo 



+ a 2 2- 



duo 
de Br 



duo 



with 



u 2 



OL\2 = 



«22 = 



jj2,Br c Br 

U \\ e \\ € ± 

U 2 (ef f 

u 2 



r,, + 7 - 



r,, + 7 - 



dT 



d'j 
dT 



-Br 



2 ' 



G?7 / ' 



U 2 (e Br fJ d l ' 

dr. r 



u±<f 



U 2 {e Br f d 7 ( e f) 2 ' 



(32) 
(33) 

(34) 

(35) 
(36) 
(37) 
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and 



3 sinh 



dlj 

g?7 




dT 



_L 



(7-1) 7 (7-1) 

3 (2 + 7) sinh -1 



1- 



7)' 



(1-7) 



g?7 



= < 



(2 + 7) sec" 1 v /7 



+ 



for < 7 < 1 

for 7 > 1 
for < 7 < 1 

for 7 > 1 



(38) 



(39) 



(7-I) 5 (7-1)' 

Now we turn to the Bruggeman equations (15) and (16). Their derivatives with respect to 
uj may be written as 

de Br de Br 
Pu-r- + Pi2-r^ + fa = , 



21- 



du 
de Br 



+ P: 



22- 



dw 
dw 



+ P23 = , 



with 



P11 
P12 

013 

(321 
022 



P. 



2.3 



an (e a - ef ) {e b - ef ) + £>„ (2ef 
ai2 (e a - ef ) (e 6 - ef ) , 

«2i (e° - ef) (e 6 - ef ) , 
« 22 (e a - ef ) (e b - ef ) + D ± (2e* 



e — e 



de b 
duj 
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The derivatives of ef and ef therefore finally emerge as 

^ 6 f /9l2/923 — P22P13 



duo 
def 



PllP22 ~ P\2$2\ 
$2\P\3 — PllP23 



(40) 
(41) 

(42) 
(43) 

(44) 

(45) 
(46) 

(47) 

(48) 
(49) 



du P11P22 ~ P12P21 

To summarize, given a uniaxial dielectric HCM with permittivity dyadic e Br estimated 
using the Bruggeman homogenization formalism, the group velocity (18) may be computed 
using the expression (26), with (27) and (28), wherein the derivatives of ef and ef are 
provided by (48) and (49). 
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4 Numerical studies 



Without loss of generality, let us choose the axis of rotational symmetry of the component 
spheroids to lie along the x axis, i.e., c — x. We consider wavevectors lying in the xy plane, 
oriented at an angle 9 to the x axis. That is, we take 

k = k{cos9,sm9,0} . (50) 

Thus, the magnitudes k = k or = |&^, r | and k = k ex = \k ex \ of the ordinary and extraordinary 
wavevectors arise from (24) as [15] 



k r = uj\//j, ef r , (51) 

(52) 



kex — ^ 



\ 



ejf r cos 2 9 + ef r sin 2 9 



Let us explore numerically the enhancement in group velocity that can arise through 
homogenization, paying special attention to directional effects induced by the shape of the 
component spheroidal particles. In particular, we choose the component material phase a 
to have a relatively high permittivity e a and a relatively small frequency-dispersion term 
de a /duj, compared with the component material phase b. As representative constitutive 
parameter values, we set: e a = 30e , (de a /duj)\ U)=ulg = 6e /u; o , e b = 1.2e„ and [de b /duo) \ UJ=U1 = 
12e /u o , where e is the permittivity of free space. 

In Figure 1, the Bruggeman estimates of the HCM permittivity parameters e^ r and e^ r 
are plotted as functions of volume fraction f a , for the range of values of p = U\\/U± shown 
in Table I. Clearly, — > e b as f a — > and e^ r y — > e a as f a — > I. We see that e^ r becomes 
an increasingly nonlinear function of f a as p decreases, whereas r becomes an increasingly 
nonlinear function of f a as p increases. 

In Figure 2, the magnitude of the group velocity Vg T = \v_g r \ of a wavepacket in the chosen 
HCM is plotted against volume fraction. The group velocities are calculated with k = k ex for 
9 = 0°, 30°, 60° and 90°. The corresponding graphs for 180° — 9 are the same as those for 9. 
Since the ordinary wavevector k^. = k ex at 9 = 0°, the ordinary group velocities for any 9 are 
identical to those provided in Figure 2(a) wherein the results for 9 = 0° are presented. The 
group velocity magnitudes for the component material phases a and b are v a g = 0.166c and 
v b = 0.152c, respectively (as is confirmed in Figure 2 by the group velocity values at f a — 1 
and f a = 0, respectively), where c = l/y/e p . Hence, for this particular homogenization 
example, group- velocity enhancement arises when Vg T > max [vg, v b g } — 0.166c. The group- 
velocity-enhancement region is identified by shading in Figure 2. 

It may be discerned from Figure 2(a) that group-velocity enhancement occurs over an 
increasingly large range of f a values as p decreases. Furthermore, the degree of enhancement 
at p = 20 is much smaller than it is at p = 0.05. 

As 9 increases, the range of f a values at which group- velocity enhancement occurs pro- 
gressively decreases for small values of p. In fact, at 9 — 60° there is no longer any en- 
hancement in group velocity for p = 0.05. At 9 = 90°, the group- velocity enhancement 
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characteristics at low and high values of p are approximately the reverse of their respective 
characteristics at 9 = 0°. That is, group- velocity enhancement occurs over a wide range of 
f a values for high values of p at 9 = 90°, but there is no enhancement in group velocity at 
low values of p. 

Clearly therefore, enhancement of group velocity is maximum in a direction parallel to the 
longest semi-axis of the spheroidal particles, which can be prolate (p < 1) or oblate (p > 1). 
For spherical particles (p = 1), group- velocity enhancement is direction-independent, and 
we recover the results of the predecessor study [9]. 

5 Concluding remarks 

The enhancement in group velocity brought about by homogenization is sensitively depen- 
dent upon directional properties. Both the shape of the component spheroidal particles, 
and their orientation relative to the direction of propagation, strongly influence the group- 
velocity enhancement. 

The homogenization scenario presented here deals with the conceptualization of a uniaxial 
HCM as arising from identically oriented spheroidal particles of isotropic component material 
phases. The homogenization of two uniaxial dielectric component phases distributed as 
spherical particles is mathematically equivalent, provided that the distinguished axes of the 
component material phases have the same orientation [16]. 

Acknowledgement: We thank Prof. G.W. Milton for drawing our attention to group- 
velocity enhancement in isotropic dielectric composite mediums with HCM permittivity es- 
timated using the Maxwell Garnett formalism, as described in Refs. [7] and [8] . 
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Table 1. Key for the values of p — U\\/U± used in Figures 1 and 2. 
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Figure 1: The HCM relative permittivity parameters e / e (above) and ef r /e (below) 
plotted against volume fraction f a . Permittivities of component material phases: e a = 30e 
and €b = 1.2e . A key for the p = U\\/U± values is given in Table 1. 
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0.2 0.4 0.6 0.8 1 

f a 

Figure 2: The magnitude of the HCM group velocity v® r = \v_g T \, as estimated using the 
Bruggeman formalism, plotted against volume fraction f a . The group velocity is normal- 
ized with respect to c = l/*/e // - Constitutive parameters of component material phases: 
e a = 30e , (de a /dcu)\ uJ=(Vo = 6e /u o , e b = 1.2e and (de b /du) \ u)=Uo = 12e /u o . A key for 
the p = U\\/U± values is given in Table 1. Shading indicates the region of group-velocity 
enhancement, (a) Extraordinary wavevector angle 8 = 0°. 
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0.2 0.4 0.6 0.8 1 

f a 

Figure 2: (c) Extraordinary wavevector angle 9 = 60°. 
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Figure 2: (d) Extraordinary wavevector angle 9 = 90°. 



